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1 Introduction. 



Parton distribution functions as gained from high-energy scattering reactions provide our main 
information about the non-perturbative structure of hadrons. Up to now the wealth of this 
knowledge was obtained from deep inelastic scattering (DIS) of leptons off a hadron target which 
allows to access the ordinary forward structure functions. Recently much excitement was generated 
by new objects which could provide a new insight into the underlying dynamics of non-perturbative 
QCD in hard processes — the non-forward parton distributions [p], @, 0- Several processes were 
proposed where the latter could be measured, e.g. deeply virtual Compton scattering, diffractive 
production of vector meson etc. It is unnecessary to emphasize that any precise analysis of 
such data will require accurate predictions from strong interaction theory for the corresponding 
reactions, especially for the radiative corrections to the tree level amplitudes. We solved the latter 
problem in a series of our recent papers (see |3J] and references cited therein) where all necessary 
ingredients for a next-to-leading order (NLO) analysis were derived. In this contribution we 
conclude by presenting a general method for the solution of the generalized evolution equations for 
the non-forward parton distributions in two-loop approximation in the flavour singlet channel. The 
technique we use, i.e. the reconstruction of a distribution function^], 0(x,(), from its conformal 
moments, Oj, via an expansion in the series of the orthogonal Jacobi polynomials ||, Pj 01 '^, 

0(x,Q<xw(x\a,p) £ P^ ) (2x-l)Y / E jk (C)O k , (1) 

j k 

is close in spirit to the one developed about two decades ago by Parisi and Sourlas || (see also 

Ref. j7|) for the usual DIS. The generalization to the singlet sector in DIS was done in Ref. |8j. 

Quite recently it was independently rediscovered by us in Refs. [|, Q in relation to the problem 

at hand. This approach is the most suitable for our purposes because the standard methods of 

contour integration via the inverse Mellin transformation are very hard to implement in practice. 

On top of this for particular values for the parameters of the Jacobi polynomials the latter become 

eigenf unctions of the leading order non-forward evolution equation which favours the choice we 

have made. Since for the time being analytical expressions in NLO exist only for the complete set 

of non-forward singlet anomalous dimensions [ |TU| , |TT| but not for the kernels themselves it is the 

only possible method to go beyond the one-loop approximation. Though, using the generating 

function for the diagonal part of the non-forward evolution kernels derived and tested at leading 

order in Ref. |T(J (the non-diagonal parts are known analytically) one can extend this procedure 

2 The functions 0{x, £) used here vanishes outside the region x £ [0, 1]. However, below we deal with non-forward 
parton distributions with the support — 1 + C < x < 1 which parametrize the matrix element of a light-ray operator. 
The results for the former case can be found in Appendix |^. 
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to two-loop approximation and find in principle NLO exclusive kernels, V(t, t') with \t\, \t'\ < 1, 
analytically. After continuation to the entire plane {\t\, \t'\ < 00} might allow a direct numerical 
integration of the two-loop evolution equations — an alternative method used in the analysis of 
the forward structure functions. 

Note, however, that with the former approach at hand, our numerical task is more difficult 
contrary to method used in usual forward scattering as we cannot restrict ourselves to the first 
iV max ~ 10 polynomials in the series (P because in DIS the shape of the curve is roughly fixed by 
the weight-function, w(x\a,(3) = x a x@ (with x = 1 — x), provided we have made an appropriate 
choice for the indices of the Jacobi polynomials, namely, with (3 given by Regge theory predictions 
and a driven by quark counting rules near the phase space boundary^. The series of polynomials 
leads only to small perturbations around this x a or -behaviour. For the case at hand the situation 
is different since the shape of the distribution is obtained from the total series over all oscillating 
polynomials. In practice, due to rather rapid convergence of the series the number of terms to be 
taken into account is rather large but still treatable. 

In the present paper we pursue the goal of developing a machinery for the solution of the two- 
loop renormalization group equations for the non-forward parton distributions and of studying 
numerically the Q 2 -dependence of the latter. Our presentation is organized as follows. In the 
next section we give our definitions and conventions used throughout the paper. The third section 
is devoted to the description of our formalism for evolution of generalized distribution functions 
together with a set of explicit formulae required in the numerical analysis. Since there are no 
experimental data yet for the non-forward parton distributions we have to rely on estimations of 
the form of the (x, £)-dependence at low normalization scale from non-perturbative approaches 
to QCD. In section 4 we give our results for models of the singlet non-forward parton functions. 
We consider CTEQ-based functions which are deduced making plausible assumptions about the 
behaviour of the so-called double distributions @] in different regions of phase space. The final 
section is left for the conclusions. To make the presentation self-consistent as much as possible 
we add three appendices. The first one contains the analytically continued forward anomalous 
dimension matrix responsible for the evolution of the multiplicatively renormalizable conformal 
moments. The second one presents the non-diagonal elements of the anomalous dimension matrix 
of the tree-level conformal operators which define the corrections to the eigenfunctions of the 
NLO evolution kernels. The last appendix contains the formulae required for reconstruction of 
the non-forward distributions with the support < x < 1. 

3 Let us remind, however, that there exists another optimal choice of a and (3 which leads to the same high 
precision reconstruction of the structure function Obviously, there is no dependence on the particular values 
of these parameters provided a large enough number of terms, N max , is taken into account. 
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2 Conventions. 



Since we are interested in the study of the parity even flavour singlet evolution equations we face 
as usual the mixing problem between quarks and gluons. To treat it in a compact way let us 
introduce a two-dimensional vector of singlet non-forward parton distributions composed from 
quark and gluon functions (( = 1 — (): 

° {xX)= {f(x()) With -C<z< L ( 2 ) 
The latter are defined as the light-cone Fourier transforms^ 

{h'\ Q 0{K U K 2 )\h) = 2 f\dxe- iK ' x - iK ^- x) Q{x,0, (3) 

(h'fOiKx, K 2 )\h) = J _dx e- iKlX - lK2 ^-' x) g{xX), (4) 

of the light-ray quark and gluon string operators {v + = v^n^) 

k 2 ) = ^(« 2 7i)7 + V'(Kin) - ^(Kira)7 + V>(K27i), g O{h^ k 2 ) = G +jU (/«2n)G M+ (/€in), (5) 

where, for brevity, we omit a path-ordered link factor which ensures gauge invariance. Here we 
accept the conventions advocated by Radyushkin ||, namely, x is the momentum fraction of an 
outgoing parton w.r.t. the incoming hadron momentum p, k + = xp + , and A + = p + — p' + = (p+, 
where p' is an outgoing hadron momentum. These quantities are related to the variables t and 1] 
used by the authors [| according to rj = ^z^, t = ^f^S where k + = tP + , A + = r]P + and the 
averaged momentum is introduced as follows P = p + p'. An advantage of the first conventions is 
that the variable x acquires a simple partonic interpretation in contrast to t. However, the range 
of the latter does not depend on the longitudinal asymmetry parameter 77 as compared to x. 

Due to charge conjugation properties of the non-local operators (|5|) it immediately follows 
that in the (x, C)-plane the singlet quark distribution is anti-symmetric, Q(x,() = —Q(C ~ x ->(), 
while the gluon one is symmetric, G(x,Q = Q{Q — w.r.t. the line x — |. Up to an obvious 
redefinition the functions Q and Q coincide with the ones introduced in Refs. [|I|. [| 

q(*= C = TX-l = (1 + V)q& 17), (6) 

and similar for gluons, where the variables t and rj vary within the limits — 1<£<1,0<7/<1. 
Obviously, these singlet quark and gluon distributions are, respectively, odd and even functions 
oft. 



4 Below, we will repeatedly omit the integration limits. The latter can be easily restored making use of known 
support properties of the distributions or/and evolution kernels. 



3 



The original functions, Q and Q, introduced above can be decomposed into the following quark, 
anti-quark, ^O, and gluon, g O, non-forward distributions defined in the range x G [0, 1] via 

Q(x,C) = 1{[ q O(xX) + Q O(xX)]0(x)6(x) 

- ~Qo((-x,() + Q O((-x,()}0((-x)e(x + ()}, (7) 

g(x,c) = l{ G o(x,oe(x)6(x) + G o((-xX)e((-x)6(x + o}. (8) 

However, our consequent discussion deals entirely with the non- forward distributions Q and Q. 

3 The method. 

Due to ultraviolet divergences of perturbative corrections for a product of operators separated by 
a light-like distance as in Eq. (|5|) the distributions acquire a dependence on a normalization point, 
H 2 , governed by a renormalization group equation, the so-called generalized Efremov-Radyushkin- 
Brodsky-Lepage (ER-BL) evolution equation 

/x 2 ^0(x, = J dx'K (x, x', C |a.) 0(x\ C), (9) 

where the purely perturbative kernel K is a 2 x 2 matrix given by a series in the coupling 
constant, a s . The LO non-forward light-cone fraction kernels^ were evaluated many years ago in 
Ref. T2] for even and odd parity and chirality while corresponding light-cone position counterparts 



were addressed in Refs. R [H|-[[r7]]. The evaluation of the two-loop corrections, however, for the 
Gegenbauer moments of the kernel K 

J dx Cf A) (2^ - AB K x>, C I a.) = -\ g AB ljk (a s )C^ B) (2^ - , (10) 

- where the numerical values of the indices v{A, B) depends on the channel under consideration 



has been addressed by us in Refs. WA P^Q , [TlH • In one- loop approximation the above kernel 



is diagonal 7-^ = 7? Sjk in this basis while beyond LO non-diagonal elements, 7^.°, appear 
7^^ = 7^. + 7jl D - The diagonal anomalous dimensions coincide up to pre-factors given in Eq. ( |T8| ) 
with the anomalous dimensions of local operators without total derivatives which appear in the 
operator product expansion (OPE) for DIS. The formalism we have developed there allowed us 
to find all entries, tJjP, in closed analytical form in NLO and provided a simple diagonalization of 
the evolution equation @. The following discussion concerns the solution of Eq. (0) within the 
formalism we have sketched above. 



3 For recent independent calculations of these kernels related to the problem at hand see the reviews |il 
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Figure 1: Two-loop diagrams giving rise to a no n- vanishing contribution to the W-functions of 
the generalized ER-BL evolution kernel QTTj) . 

Let us add a remark about the general structure of the all-order evolution kernel riV (-, — ] = 
K (x = j^j, x' = C = T^rj) • ^ can ^ e deduced in a straightforward manner from the support 
properties and charge conjugation symmetry of the light-ray operators and reads, e.g. for the 
QQ-channel |]§, 

v(t, t>) = e° u (t -t',t- i)v(t, t>) + e° n (-t - 1>, -t - 1) w(-t, + { *Z } » ( n ) 

where we have introduced the step-function e^t.f) = ^ {^(t)6>(— - Q and 

V, W are analytic functions of their arguments. Here the W-part appears first in NLO and is 
generated by diagrams depicted in Fig. |l[ These terms give rise to mixing of partons from different 
regions of phase space conserved by the leading order evolution and it forces us to consider the 
functions ([7],^) which are defined on the entire interval of x G [— (, 1]. 



3.1 General formalism. 

In order to solve the evolution equation (Q) one has to find its eigenvalues and eigenfunctions. 



This problem has been exhaustively treated by us in general form in Ref. |TT| where it was shown 
that its solution can be written in terms of the partial conformal wave expansion 



0(x, C, Q 2 ) = E 4>j C |«.(Q 2 )) Oj(C, Q 2 



3=0 



with the partial conformal waves matrix 



i (x, ( |« S (Q 2 )) = E C fc - J fc (x, B kj , 

k=j 



(12) 



(13) 
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being the eigenstate of the all-orders equation (Q). The B- matrix defines the corrections to the 
eigenf unctions which at tree-level diagonalize the leading order ER-BL equations. The latter readf] 

(w(x\u) = (xx) v ~ l/2 ) 



4>i (x, = 7 



1 I ?^|)^ V 2 ? ~ l ) u 



(14) 



The multiplicatively renormalizable moments, Oj((,Q 2 ), evolve as follows 

d,(C,g 2 ) = £ j (a s (Q 2 ),a s (Q 2 )) O (15) 
with an evolution operator defined by the equation 

^(a.CtfWQg)) = Tex p{-^/ Q T V 7 ? (as(r)) } (16) 
where the operator T orders the matrices of the diagonal anomalous dimensions 



IQl^D qg d 

7? = I I (17) 



ij 



J 7 



along the integration path. They are given as an expansion in the coupling constant by "y • (a s ) = 
'Y*; ' 1 + ff^ 'v!; + . . . — The elements of this matrix are related to the forward anomalous 

ZTY ' J \ ZTV J ' J 

dimensions 7^ w via the relations 



QQrS> _ QQ^iw qg d _ 6 qg fw god = 3_ gq > fw ggd = GG fw 
n H ' ij j H ' <i q ij ' 'i ij 



The pre-factors 6/j, j/6 in the off-diagonal matrix elements come from the conventional definition 
of the Gegenbauer polynomials. 

The normalization condition in Eqs. (|T^,[T^) is chosen so that there are no radiative corrections 
at an input scale Qq. Apart from the minimization of the higher order corrections the advantage 
of this choice lies in the fact that multiplicatively renormalizable Oj(CQg) are defined at Qo by 
forming ordinary Gegenbauer moments with non-forward distributions via 



J ((,Ql)= jdxC 3 {x,QO{xX,Ql), withC^C) 



' C J C/ (2f - l) 

v <?- x c)_ x (2f - 1 



(19) 



6 Let us emphasize that here the Gegenbauer polynomials should be understood as mathematical distributions 
(see Ref. Q) according to the relation Q ±w(f\v)Cj (2f — l) oc J* dy(yyy +u ~ 1 ^ 2 S^ ((y — x) in order to be able 
to restore the correct support properties of the non- forward distributions. 
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As before [|J the restoration of the support properties of the distributions is achieved via an 
expansion of the latter in a series with respect to the complete set of orthogonal polynomials, 
•p (Qp) (t)> on the interval -1 < t < 1 



0(x,CQ 2 



2 E^f^Wc,^ 



s .7=0 



2-C 







/ w(t\ a p) 'p( a p) ( f\ 



nj(a' p ) 1 5 W 



(20) 



(21) 



with w(t\a p ) and nj(a p ) being weight and normalization factors, respectively (see, e.g. [[H|]). The 
matrix of moments, Alj(C, Q 2 ), is given by the sum 



Mj(C,Q 2 ) = £s jfc (C)O fc (C,Q 2 ), where 



k=0 
3 



(22) 
(23) 



k=0 



and the upper limit in Eq. (|2~2"1) will come from the constraint ^-functions in the expansion coeffi- 
cients 







E jk (C) = 

with elements given by the integral (where 6j k = {1, if j > k; 0, if j < k}) 



(24) 



jk T(v)(u + k) f i A+ n +2 . k+u , d k (ap) ( (t 



, dt (l- t z r +u ~2— v 

(20 k Ta)T(k + u+ 1 1 )J-i [i t] dt k ^ 



2-C, 



(25) 



Taking Jacobi polynomials, Vj ap \t) = Pj a ' p> (t), we obtain the rather complicated result 



E J jk {v-a,(3\Q = 



^ Tjy) r(j + q + l) T{k + v+\) 
T{\)T{j + a + (3+l) T(2k + 2u) 

T(j + k + 1 + a + /3 + 1) 



.(2_ C )-*-V\ 



2fc + 2u + 1 



(26) 



T(£ + l)T(j — k — £+ l)T(k + £ + a + 1) 

which significantly simplifies for the particular values of the indices a = (3 = // — ~ when Jacobi 
polynomials degenerate into Gegenbauer ones. Therefore, taking the latter as a set of expansion 
functions for our purposes, i.e. v\ ap \t) = Cj(t) we have 



T(u) 



(-i)^r „ + ^ 



r(/i) r (z/ + fc) r ( 1 



/ _i^e / 1 4. it* 
v I i/ + A; + 1 



(27) 



(2-C) 2 



The fact that only even j — /c-moments contribute is obvious since the non-forward distributions 
are even functions of rj = due to time-reversal invariance and hermiticity Jj]], — as a result 
only even powers of r\ can appear in the expansion. 

Now we are in a position to give results for the evolution operator, £, and .B-matrix in two- loop 
approximation. 

3.2 Explicit NLO solution. 

In NLO the evolution operator satisfies the equation 

dl J sm £ i {<Q 2 )MQD) = -j o {7? (0) + {*s(Q 2 ),a s (QD) , (28) 

where 

and the boundary condition £j (a s (Ql),a s (Ql)) = t with the unit matrix 1 = (qi) - Here (3q 
and (3\ are the first and second coefficient in the perturbative expansion of the QCD /3-function 
f = % A, + (^) 2 /9i + ... and read ft = %T F N f - fC A and fa = fC A N f + 2C F N f - f C 
respectively. 
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The solution of the above equation is p2|, |23 



Sj (a s (Q 2 ),a s (Ql)) = [ P+ - ^ ^£^M - P + P . P+ 1 f " ^ 



a s (Q 2 ) - a s {Ql) 1 ^ \ (<*,((% 



q 2 ) /vm'; A /« s (^)\ (77 " 7;+/3o)/A ^ (<XB(Qi)y' lfh 



2vr 77~7i+A)\ V^(Q 2 )/ j\oc 8 {Q 2 )) 
+ (+<->-), (30) 



where we have introduced projection operators 



^ = ^(7? ( °'-77l), (3D 



is -1, 

and the eigenvalues of the LO anomalous dimension matrix 

7 ± = 1 ^(0) + GC 7 D(0) ± ^Q^D(O) _ GGy D(0)^ - 4Gay D(0) QGy D(0)^ ^ 

The .B-matrix which fixes the corrections to the eigenfunctions of the NLO ER-BL kernels is 
given by 

B = l + B {1) , (33) 



S 



where is determined by the first order differential equation |TT 

d B« (a s (Q 2 ), a s (Ql)) = -J Q { [7 D(0) , B« (a s (Q 2 ), a s (Q, 



cilna s (g 2 )- 
and reads (j > fc) 



+ 



2tt 



"7 



ND(1) 



+ (+--). 



(34) 



(35) 



The two-loop forward anomalous dimension matrix in the singlet sector has been evaluated 
in Refs. |25|. The non-diagonal entries, 'y ND ( 1 ) ; of the full anomalous dimension matrix 7 = 
7° + 7 ND of the conformal operators have become available quite recently flO, IT|. For the reader's 
convenience both of the above ingredients are summarized in appendices and 0, respectively. 

Finally, the coupling constant in the corresponding order is given by the following inverse-log 
expansion 



a s (Q 2 ) 



1 + 



A In ln(Q7Afjs, 



MS , , $ hq 2 /Ay 



(36) 



Note, however, that due to the fact that the input scale, Qo, — at which the initial conditions 
considered below are defined — is very low it might be more reliable and accurate to obtain a s in 
two-loop approximation by solving the exact transcendental equation 

Q 2 4vr ft, / 4vr /V 



-A) In 



MS 



_ HI 1 

a s (Q 2 ) fa n 



47T 

'/3 a s (Q 2 



(37) 



This leads to ~ 15% discrepancy between them for Q 2 = 0.6 GeV 2 which goes down to 3% for 
Q 2 = 4 GeV 2 . However, we leave aside this source of theoretical uncertainty, as being conceptually 
irrelevant for our study it presents only an extra source of unnecessary complication. 



4 Evolution of the model distributions. 

For practical applications we have chosen the indices of the Jacobi polynomials in the simplest 
way, i.e. a — /3 — 0. In this case the former coincide with Legendre polynomials, Pj(t) = 
pj°'°\t) = Cj /2 {t) with w(t) = 1 and rij = (2j + l)/2 in Eqs. ©fn]). From the last equality 
we can read off explicit expressions for the expansion coefficients Ejf. in Eq. (|27D. Contrary to 
our previous procedure advocated in Ref. B where the evolution of a limited number of points 



9 



in the discrete {x, C}-plane has been performed exactly and interpolated afterwards by a smooth 
functions in spline approximation, presently we evolve the expansion coefficients of the non-forward 
parton distribution function in the series over orthogonal polynomials and then sum them back. 
Therefore, no smoothing recipes are required but the result depends on the truncation parameter 
iV mm . To achieve rather high reconstruction accuracy, i.e. ~ 10~ 2 — 10~ 3 , we have retained up to 
iV max = 100 polynomials in the series^] (p0|) . Even better precision, 10~ 4 , is feasible by doubling 
the number, iV max = 100, of polynomials in the expansion at the price of the much more time 
consuming procedure. Note that only the region around the crossover point x = ( is sensitive to 
JV max due to rather rapid change of the curve on the small interval of x. Taking only iV max ~ 50 
the reconstructions accuracy decreases down to 10 _1 — 10~ 2 . 

Note that due to symmetry properties of the distribution functions only those moments survive 
which acquire an operator content and evolve with anomalous dimensions deduced from OPE. This 
is in contrast to the case of distributions, "^O, with the support < x < 1 where all moments 
enter on equal grounds due to lack of charge conjugation symmetry in the parametrization of the 
hadronic matrix elements in Eqs. (H,Ei,0ll). They were considered in Ref. @] and it was crucial 
there to use the analytically continued anomalous dimensions^ (see appendix [A]) since otherwise 
the convergence of the series will break down already for iV max ~ 10. This is so in spite of the 
fact that the relative difference between the former and the ones deduced from an OPE analysis 
with a = (— iy + (see appendix [Ap left intact is negligible (< 10~ 4 ). The reason is the factorial 
growth of the expansion coefficients which weight the contribution of moments. 

Below we report on our study of the Q 2 -dependence for the distribution functions proposed in 
Ref. [^6j. We will, however, neither speculate on the physical relevance of the latter nor discuss 
their adequacy to the real world since this is a disputable issue in the lack of any experimental 
data, but we rather accept them in order to test our formalism. We have chosen for our analysis 
N f = 4 and = 246 MeV. 

J Mb 

The non-forward functions ^(x, () in Eqs. (0,^|) for the parton species A = Q, Q, G are 
defined in terms of the double distribution function A F(y, z) via the following relation || 

^(x, C, Ql) = t dy ^ dz A F(y, z, Ql) 9(1 - y - z)S(x -y-(z). (38) 
Jo Jo 

The functional dependence of A F in the y-subspace is given by the shape of the forward parton 
density while its |-dependence has to be similar to that of the distribution amplitude. This results 

7 The C++ code used in the calculation is available via ^ittp: / /www. physik.uni-regensburg.de/~nill7791 . 

8 In the above mentioned paper we have expanded the anomalous dimensions in the series = m ci rn 
with the first few terms kept in the expansion which provides a highly accurate approximation, sufficient for 
numerical studies. 
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Figure 2: Evolution of the non-forward singlet quark distributions Q(x, £). The input function 
at Q = 0.7 GeV is shown by the short-dashed line at different £'s. The full curves moving away 
from the initial function correspond to LO results for Q 2 = 10, 10 2 , 10 14 GeV 2 , respectively. The 
long-dashed lines give the NLO results for the same values of the momentum scale in the same 
order. 
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in the following model for A F(y, z, Qq) 

A F{y,z,Cfi)=A(y,Ql)Mv,z), (39) 
with the plausible profiles A 7r(y, z 



2 

Q ir(y, z) = 6- 3 (y-z), G n(y, z) = 30^(£ - z)\ (40) 

y.i yb 

for quarks and gluons, respectively. In Eq. (|39|) the function A(y, Qq) = q(y), q(y), yg(y) is 
an ordinary forward parton density measured in deep inelastic scattering at a low normalization 
point. We use the CTEQ4LQ parametrization |27[] of the parton densities at the momentum scale 
Qo = 0.7 GeV: 

xu v {x) = x(u - u)(x) = 1.315x a573 (l - x) 3 - 281 (l + 10.614a; L034 ), 
xd v (x) = x(d - d)(x) = 0.852x a573 (l - x) 4 - 060 (l + 4.852x - 693 ), 

x(u + d)(x) = 0.578x ai43 (l-x) 7 - 293 (l + 1.858x L000 ), (41) 
xg(x) = 39.873s 1 - 889 (l-x) 5 - 389 (l + 0.6 1 8x - 474 ). 

The results of the evolution are shown in Figs. ^ for singlet quark and gluon distributions, 
respectively. Due to the low input momentum scale, Qq, the change of the shape is very prominent 
already for Q 2 = 10 GeV 2 . We have plotted the curves obtained using the LO (full lines) and 
NLO (dashed lines) formulae. The difference between them is especially sizable (~ 10 — 30%) for 
quarks at small £'s and moderately large Q 2 . The limiting lines correspond to the extremely large 
momentum scales, Q 2 = 10 14 GeV 2 , when the distributions approximately reach their asymptotic 
forms (Q 2 — > oo) [p6[| : 

C 2 V" w V'c, "J ' " " ^'^ ^ ^ V 



a,(.r. o x ^ ( i - ~) (i- - i) e(c - x), £ as (*, o cx ^ fi - e(c - x). (42) 



5 Conclusions. 

In the present paper we have generalized the formalism for the solution of NLO evolution equations 
for the non-forward parton distributions developed by us previously for the non-singlet case 0, [J 
to the flavour singlet channel. It is based on the use of the Jacobi polynomials to reconstruct the 
(x, C)-dependence of the function from the multiplicatively renormalizable conformal moments 
which diagonalize the generalized ER-BL evolution equation. Due to relatively fast convergence 
of the series in orthogonal polynomials we limit ourselves to at most iV max < 100 terms. The 
resummation can still be handled numerically and the precision achieved is sufficiently high for 
any practical application. Let us note that in absence of analytic expressions for the complete 
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NLO non-forward kernels in the singlet channel (at least for the time being) there is no feasible 
alternative to the procedure described above. 

We did not analyze in full, however, how the choice for the parameters a and f3 influences the 
convergence properties of the series. One should note that a clever choice can lead to an increase 
of reconstruction accuracy by an order of magnitude for the same iV max M . 

We have studied the general pattern of evolution obeyed by the models for the non-forward 
distributions proposed recently. We have found that the difference between LO and NLO does 
not exceed 10 — 30% which suggests that the latter effects could be taken into account provided 
very high accuracy data points will be available which seems to be a very hard experimental task. 

An application of our considerations to the evolution of the parity-odd densities is straight- 
forward since the non-diagonal elements, 7 ND , of the corresponding anomalous dimensions are 



known from Ref. ]TTJ while the analytically continued two-loop forward entities [fJHJ were derived 
in JZJJ. 

We acknowledge helpful discussions with A.V. Radyushkin and would like to thank him for 
reading of the manuscript and useful comments. This work was supported by BMBF and the 
Alexander von Humboldt Foundation (A.B.). 

A Diagonal (forward) anomalous dimensions. 

Let us add few remarks concerning the NLO anomalous dimensions for DIS. It should be noted that 
beyond leading order the moments of the DGLAP splitting kernels coincide with the anomalous 
dimensions available in the literature evaluated by OPE methods only for even (odd) moments 
provided we are interested in the crossing odd (even) combination of quark-ant i- quark species. We 
give below analytically continued anomalous dimensions from even (odd) to all values of moments 



j [3^, |23[ where the above unfortunate feature is naturally withdrawn. 



The one-loop quantities are well known and read [[31] 



Q K i0) = -^(-^ + 2) + 4^(l) + , + i) 2 , + 2) + 3), (A.1) 

QG 7r (0) = -AN f T F , C +3j 't, 4 „ (A.2) 

h f 0' + l)(j + 2)(j + 3)' 1 ' 

- -c* Ho- + o + + ^^+3) - a ■ <-> 
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The anomalous dimensions at 0(a 2 s ) are given by [|25| : 
^ wNS( V) = fa - -C F C A ) I , . 4( 'i + 3) »SU + 1) - 2 3j * + + UJ + 3 



2 v [(j + iy( j + 2 y - ' ( J + i)3 (j + 2 )3 

+ 4 (^(j + 1) - ( - + 1) 1 (j + 2) ) (S 2 (j + 1) - + 1)) 
+ 16S(j + 1) + 6<S 2 (j + 1) - 7 - 2^0' + 1) + 4a- 2 " r ()/ ' ^ 



4 dV ' !./ • Q 3 (j I 2f 



9 (J + !)=»(;+ 2)3 



- ^o-+D^+D+^+i)(-f + (j + 1) 2 (j+2) ) 

43 1 151j 4 + 867j 3 + 1792j 2 + 1590j + 523 1 
24~9 (j + 1)3(j + 2)3 J 

+ ^{-^c^^^D.i.i^t^}. ,a,) 

= ^r s "v = i) - 4 w + 5 X + + 2 St 4 2 w ts7 + 1M < A - 6 > 

<Wyf « = -2C A T F N, { (-2S?CJ + 1) + 2*0 + 1) - 2SM + 1)) {J+ ( ) ^^ + 3) 
J" + 15/ + 99/ + 382/ + 963/ + 1711/ + 2273/ + 2252/ + 1488j + 480 



+ 2- 



j<J + lHj + 2)3(j + 3)3 



2C F T F N f | (2SIU + 1) - 2S 2 (j + 1) + 5 



(j + 1)(j + 2)(j + 3) 



SiU + 1) 11j 4 + 70j 3 + 159j 2 + 160j + 64 | 

(J + 1) 2 0' + l) 3 a + 2) 3 (j + 3) /' 1 ' J 



-7, 



GQ - M1) = { (-2S!(j + 1) + 10^(j + 1) - 2S 2 (j + 1)) ^ + + 1 )( J + + 4 2) 

5i (j + 1) 12j 6 + 102j 5 + 373j 4 + 740j 3 + 821 j 2 + 464 j + 96 ' 



(j + 2) 2 j(j + l) 3 (j + 2) 3 J 

j 2 + 3j + 4 
j(j + l)(j + 2) 

1 109j 9 + 1602j 8 + 10292/ + 38022j 6 + 88673j 5 + 133818j 4 + 128014/ 



- 2C A C F \ (Sl{j + 1) + S 2 (j + 1) - S' 2 (j + 1)) -rf. 
+ 



2 



9 j 2 (j + 1) 3 (j + 2) 3 (j + 3) 

1 72582j 2 + 21384 j + 2592 1 17j 4 + 68j 3 + 143j 2 + 128j + 24 
+ 9j 2 (j + 1) 3 (j + 2) 3 (j + 3) 2 "3 f(j + l) 2 (j + 2) } , 
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= |-^^o- + 1) + 1 + 1 19i4 + 114j ' 3 + 275i2 + 312i + 138 ^ 



+ C F T F N f \2 + A 



9 v ' 3 9 iO' + l) 2 (j + 2) 2 0' + 3) 

2 j6 + 16j -5 + 51j -4 + 74j -3 + 41j -2 _ gj _ 16 ^| 



+ C A < —Siij + 1) + 165i(j + 1)- 



9 iW ' AW ' J 2 (j + l) 2 (j + 2) 2 (j + 3) 2 

« 5f a. n q/ / ■ i i \ 1 457j 9 + 6855j 8 + 44428j 7 + 163542/ 
+ 85(j + 1) -S 3 (j + 1) 



9 i 2 0' + l) 3 (j + 2) 3 (i + 3)3 
1 376129j 5 + 557883j 4 + 529962j 3 + 308808j 2 + 101088j + 15552 1 
9 W+WU+WU+W J " 



(A.9) 

Here we should use the following expressions for the analytically continued functions |23| 



Si(j) = 7s + ^0' + l), (A.10) 

S 2 (j) = c(2)-v'(i + i), (A.11) 

Ssti) = C(3) + ^"(j + l), (A.12) 

S'eti) = + + (A ' 13) 



SU) - 4c(3) + .(fi-fl(,(^)-,(0) + />^K,A,4, 



where V'^^O) = d( is the poly-gamma function and 7^ = —"0(1) is the Euler-Mascheroni 
constant. 

The integral over the dilogarithm Li2(a;) can be evaluated in terms of an analytical function 



of j by using the least square fit approximation [ 23| : 

~ 0.0030 + 1.0990s - 1.5463a; 2 + 3.2860a; 3 - 3.7887a; 4 + 1.7646a; 5 . (A.15) 

1 + x 

The a = (— 1) J+1 takes the following values for "non-singlet" — charge conjugation odd, q v , and 
flavour NS, g NS , — combinations: 

a = — 1 for q v = q — q, a = 1 for g NS = [u + u) — (d + d), etc., (A. 16) 

while the singlet distributions are evolved with a = 1. 



16 



B Non-diagonal anomalous dimensions. 



The elements of the matrix 7^? derived in Ref. |TT| read 



QQ ND(1) 
Ijk 



QG ! ND(1) 



(B.l) 



(B.2) 



, QG L D (°)GG 

+ 7j ^j* 



i/jfc 7fc 5 



GO ND(1) 
Ijk 



G S, D< °' - G Sr') (a - «tf 



)(0) 



,gqD(o)qq GG G q D(0) / GG D(0) QQ D(0)^ G q 

t Ij 9jk ~ 9jk Ik + i 7j — 7fc 



'gG °(°) GG D(0)\ , goD(O) ,„ . 

7j - 7fc J 7 fc (B.3) 

- £K *£ (0) ) + °°9»} (B.4) 

Here the leading order diagonal anomalous dimensions are given by Eqs. (U) and ( [All) ; the d 
and g elements are 



GG r ND(1) 

7jk 



GG °(0) _ GG D (0) 

7j 7fe 



(2fc + 3) 



(j-A;)(j + fc + 3)' 



and 



«% fc =-c F [i + (-ir fe 



(3 + 2fc) 



x < 2Ajk + (Ajk — + 2) + 



U-k)(j + k + 3) 

(j-k){j + k + 3) 



(k + l)(k + 2) 



GQ 



'9jh 



-C F 



j-k 



(3 + 2k) 



GG, 



9jk 



-c A [i + (-iy~ k ] 



fj-2,k 



x !2A jk + (A jk -iP(j + 2)+iP(l)) 



6(fc + l)(fc + 2)' 
(3 + 2A;) 
(j-k)(j + k + 3) 

T(j+4)r(A;) 



(B.5) 

(B.6) 
(B.7) 
(B.8) 
(B.9) 



row + 4) 



+ 2{j-k){j + k + 3) 



T(k) 



r(fc + 4) 



respectively. We have introduced here the matrix A whose elements are defined by 



.4 



1> 



'j + k + ^ 



'3 ~ k" 



+ 2i/>(j-k)-il>(j + 2)-,l>{l). 



(B.10) 



C NFPD with support < x < 1. 

Here we give a list of formulae for the expansion of non-forward parton distributions with the 
support < x < 1. The only difference which arises w.r.t. the results given in the body of the 
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paper is the argument of the Jacobi polynomials, PL- '{2x — 1), in the expansion 



0(x,(,Q 2 ) = J2'P 3 (x)M J (C,Q 2 ) 



(C.l) 



3=0 



with elements V (x) = P j "'^ ( 2x _ X ) where n j( a i P) = ^jla+p+i)mt+ali+i) • The Jacobi 



moments are given by Eq. (f22|) provided we will substitute the elements of the coefficient matrix 



The results for other classic orthogonal polynomials can be immediately derived from this equati 
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